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TYPE  19  FLOOD  INSURANCE  STUDY:  TSUNAMI  PREDICTIONS 
FOR  SOUTHERN  CALIFORNIA 


PART  I :  INTRODUCTION 

Background 

1.  Of  all  water  waves  that  occur  in  nature,  one  of  the  most  de¬ 
structive  is  the  tsunami.  The  term  "tsunami,"  originating  from  the 
Japanese  words  "tsu"  (harbor)  and  "nami"  (wave),  is  used  to  describe  sea 
waves  of  seismic  origin.  Tectonic  earthquakes ,  i.e.,  earthquakes  that 
cause  a  deformation  of  the  seabed,  appear  to  be  the  principal  seismic 
mechanism  responsible  for  the  generation  of  tsunamis.  Coastal  and  sub¬ 
marine  landslides  and  volcanic  eruptions  also  have  triggered  tsunamis. 

2.  Tsunamis  are  principally  generated  by  undersea  earthquakes  of 
magnitudes  greater  than  6.5  on  the  Richter  scale  with  focal  depths  less 
than  30  miles.*  They  are  very  long-period  waves  (5  min  to  several  hours) 
of  low  height  (a  few  feet  or  less)  when  traversing  water  of  oceanic  depth. 
Consequently,  they  are  not  discernible  in  the  deep  ocean  and  go  unnoticed 
by  ships.  Tsunamis  travel  at  the  shallow-water  wave  celerity  equal  to 
the  square  root  of  acceleration  due  to  gravity  times  water  depth  even  in 
the  deepest  oceans  because  of  their  very  long  wavelengths.  This  speed 

of  propagation  can  be  in  excess  of  500  mph  in  the  deep  ocean. 

3.  When  tsunamis  approach  a  coastal  region  where  the  water  depth 
decreases  rapidly,  wave  refraction,  shoaling,  and  bay  or  harbor  reso¬ 
nance  may  result  in  significantly  increased  wave  heights.  The  great 
periods  and  wavelengths  of  tsunamis  preclude  their  dissipating  energy 
as  a  breaking  surf;  instead,  they  are  apt  to  appear  as  rapidly  rising 
water  levels  efad  only  occasionally  as  bores. 

U.  Over  500  tsunamis  have  been  reported  within  recorded  history. 
Virtually  all  of  these  tsunamis  have  occurred  in  the  Pacific  Basin.  This 


*  Multiply  miles  by  1.6093UU  to  convert  to  kilometres. 
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is  because  most  tsunamis  are  associated  with  earthquakes,  and  most 
seismic  activity  beneath  the  oceans  is  concentrated  in  the  narrow  fault 
zones  adjacent  to  the  great  oceanic  trench  systems  which  are  predomi¬ 
nantly  confined  to  the  Pacific  Ocean. 

5.  The  loss  of  life  and  destruction  of  property  due  to  tsunamis 
have  been  immense.  The  Great  Hoei  Tokaido-Nankaido  tsunami  of  Japan 
killed  30,000  people  in  1707.  In  1868,  the  Great  Peru  tsunami  caused 
25,000  deaths  and  carried  the  frigate  U.S.S.  Waterlee  1,300  ft  inland. 

The  Great  Meiji  Sanriku  tsunami  of  1896  killed  27,122  persons  in  Japan 
and  washed  away  over  10,000  houses. 

6.  In  recent  times,  three  tsunamis  have  caused  major  destruction 
in  areas  of  the  United  States.  The  Great  Aleutian  tsunami  of  1946  killed 
173  persons  in  Hawaii,  where  heights  as  great  as  55  ft  were  recorded. 

The  i960  Chilean  tsunami  killed  330  people  in  Chile,  6l  in  Hawaii,  and 
199  in  distant  Japan.  The  most  recent  major  tsunami  to  affect  the  United 
States,  the  1964  Alaskan  tsunami,  killed  107  people  in  Alaska,  4  in 
Oregon,  and  11  in  Crescent  City,  California,  and  caused  over  100  million 
dollars  in  damage  on  the  west  coast  of  North  America. 

Purpose  of  Study 

7.  The  purpose  of  this  study  was  to  establish  100-  and  500-yr 
tsunami  elevations  in  southern  California  produced  by  distantly  generated 
tsunamis.  This  study  is  an  update  of  a  previous  study  that  predicted 
100-  and  500-yr  tsunami  elevations  in  southern  California  (Houston  and 
Garcia  1974).  Whereas  the  previous  study  by  Houston  and  Garcia  (1974) 
used  simple  analytical  solutions  to  propagate  tsunamis  over  the  nearshore 
region  and  to  statistically  combine  tsunamic  and  astronomical  tides, 
this  study  used  complex  numerical  solutions  for  nearshore  propagation  and 
for  superposition  of  tsunamis  and  tides.  Two  other  reports  (Garcia  and 
Houston  1975  and  Houston  and  Garcia  1978)  established  tsunami  elevations 
for  the  west  coast  of  the  continental  United  States  excluding  Southern 
California.  The  100-  and  500-yr  tsunami  elevations  are  required  by  the 
Federal  Insurance  Administration  (FIA)  of  the  Federal  Emergency  Manage¬ 
ment  Agency  (FEMA)  for  use  in  flood  insurance  rate  calculations. 
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PART  II: 


NUMERICAL  MODELS 


Background 

8.  Unlike  other  areas  of  the  Pacific  Ocean  such  as  the  Hawaiian 
Islands,  southern  California  lacks  sufficient  data  to  allow  reasonable 
tsunami  elevation  predictions  based  upon  local  historical  records  of 
tsunami  activity.  Virtually  all  of  southern  California  is  completely 
without  data  of  tsunami  occurrence,  even  data  for  the  prominent  i960 
and  196U  tsunamis.  Only  a  handful  of  locations  have  historical  data  for 
tsunamis  other  than  the  i960  and  196**  tsunamis.  However,  the  FIA  re¬ 
quires  information  on  tsunami  elevations  for  all  of  southern  California, 
even  for  the  many  locations  that  have  no  known  historical  data  of  tsunami 
activity  and  for  coastal  areas  that  are  currently  not  developed  (since 
these  areas  may  be  developed  in  the  future). 

9.  The  lack  of  historical  data  of  tsunami  activity  in  southern 
California  necessitates  the  use  of  numerical  models  to  predict  tsunami 
elevations.  The  Aleut ian-Alaskan  area  and  the  west  coast  of  South 
America  were  found  by  Houston  and  Garcia  (197*+)  to  be  the  tsunamigenic 
regions  of  concern  to  southern  California.  Both  regions  have  sufficient 
data  on  the  generation  of  major  tsunamis  to  allow  a  statistical  investi¬ 
gation  of  tsunami  generation.  A  numerical  model  employing  a  fairly 
coarse  grid  to  cover  a  large  section  of  the  Pacific  Ocean  was  used  to 
generate  representative  tsunamis  and  propagate  them  across  the  deep 
ocean.  A  second  numerical  model  employing  a  fine  grid  was  used  to  propa¬ 
gate  tsunamis  from  the  deep  ocean  over  the  continental  slope  and  shelf 

to  shore.  Previous  studies  by  Houston  and  Garcia  (197*0  and  Garcia  and 
Houston  (1975)  used  analytical  methods  to  solve  one-dimensional  linear 
equations  that  described  the  transformation  of  tsunamis  as  they  propagated 
from  deep  water  to  shore.  Houston  and  Garcia  (1978)  used  a  two- 
dimensional  numerical  model.  In  this  study  a  two-dimensional  model  is 
used  for  nearshore  propagation  that  is  similar  to  the  model  used  by 
Houston  and  Garcia  (1978)  but  allows  a  numerical  grid  with  a  varying 
spatial  grid  cell  size. 
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10.  There  has  been  disagreement  among  tsunami  researchers  concern¬ 
ing  the  equations  that  govern  tsunami  generation,  deep-ocean  propagation, 
and  nearshore  propagation.  The  significance  of  frequency  dispersion  and 
nonlinear  advection  on  tsunami  propagation  has  been  determined  only  in 
recent  years.  Using  the  nonlinear  and  dispersive  Korteweg-deVries  equa¬ 
tion,  Hammack  and  Segur  (1978)  concluded  that  "if  the  length  of  the 
initial  wave  is  approximately  100  miles,  the  lead  wave  is  described  by 

a  linear  nondispersive  model  from  source  region  until  shoaling  occurs 
near  the  shoreline."  This  conclusion  also  applies  to  leading  waves  that 
have  reflected  off  land  areas  and  traveled  to  a  distant  location.  Thus 
two  or  three  of  the  initial  waves  may  be  governed  by  linear  nondispersive 
equations.  Similarly,  in  studying  tsunami  propagation  from  the  deep 
ocean  to  the  nearshore  region.  Goring  (1978)  concluded  that  "...  because 
of  the  small  relative  height  of  tsunamis  and  their  large  lengths  relative 
to  the  lengths  of  the  continental  slope,  the  propagation  of  tsunamis  from 
the  deep  ocean  to  the  continental  shelf-break  and  for  some  distance  onto 
the  shelf  will  be  predicted  as  well  by  the  linear  nondispersive  theory 
as  by  the  nonlinear  theories."  Furthermore,  Tuck  (1979)  found  that 
"...  linear  long-wave  equations  are  adequate  to  describe  most  of  the 
tsunami  generation,  propagation  and  reception  processes."  Also,  numeri¬ 
cal  studies  by  Hwang  et  al.  (1972)  have  indicated  that  linear  nondisper¬ 
sive  equations  govern  tsunami  propagation  over  the  deep  ocean.  Numeri¬ 
cal  studies  by  Houston  (1978)  have  shown  that  linear  nondispersive  equa¬ 
tions  govern  tsunami  generation,  propagation  over  the  deep  ocean,  and 
interaction  with  the  Hawaiian  Islands. 

11.  In  this  study,  only  large  tsunamis  (initial  wavelengths 
greater  than  100  miles)  are  considered;  therefore,  linear  nondispersive 
equations  are  used  for  tsunami  generation  and  deep-ocean  propagation. 

The  nearshore  numerical  model  includes  bottom  friction  and  nonlinear  ad¬ 
vection  terms.  However,  it  was  found  by  numerical  experiments  that  these 
terms  have  no  significant  effect  on  tsunamis  propagated  from  the  deep 
ocean  to  the  shoreline  in  the  southern  California  region. 

12.  In  this  report,  only  tsunamis  of  distant  origin  (at  least  two 
or  three  tsunami  wavelengths  away  from  the  west  coast )  are  considered  in 


the  analysis.  Hammack  (1972)  has  shown  that  near  the  generation  area  of 
a  tsunami,  details  of  ground  motion  during  the  earthquake  and  details  of 
the  permanent  deformation  of  the  sea  floor  influence  the  form  of  the  re¬ 
sulting  tsunami.  Very  little  is  known  about  the  actual  time-dependent 
ground  motion  during  earthquakes  generating  major  tsunamis,  and  small- 
scale  details  of  the  permanent  deformation  of  the  sea  floor  following 
earthquakes  cannot  be  predicted  in  advance.  Thus,  accurate  predictions 
of  the  properties  of  xceally  generated  tsunamis  are  not  possible  at  this 
time.  However,  Hammack  (1972)  has  shown  that  the  time-dependent  ground 
motions  and  small-scale  details  of  the  permanent  ground  deformation 
produce  waves  which  are  not  significant  far  from  the  source  region. 

Thus,  distantly  generated  tsunamis  can  be  studied  only  knowing  major 
features  of  the  permanent  ground  deformation. 

13.  The  probability  is  not  considered  great  that  an  extremely 
destructive,  locally  generated  tsunami  will  occur  in  southern  California. 
Tsunamis  are  generally  produced  by  earthquakes  having  fault  movements 
that  exhibit  a  pronounced  "dip-slip,"  or  vertical  component  of  motion. 
"Strike-slip,"  or  horizontal  displacement,  fault  movements  are  ineffi¬ 
cient  generators  of  tsunamis.  Faults  in  southern  California  character¬ 
istically  exhibit  strike-slip  motion  since  the  Pacific  block  of  the 
earth's  crust  is  moving  horizontally  relative  to  the  North  American  block. 

14.  There  have  been  reports  of  significant  locally  generated 
tsunamis  in  southern  California.  For  example,  a  recent  publication  of 
the  California  Division  of  Mines  and  Geology  (Weber  and  Kiessling  1978) 
mentions  that  Wood  and  Heck  (1966)  reported  that  runup  heights  of  a 
tsunami  generated  by  the  1812  Santa  Barbara  earthquake  reached  50  ft*  at 
Gaviota,  30-35  ft  at  Santa  Barbara,  and  15  ft  or  more  at  Ventura  in 
southern  California.  However,  an  exhaustive  study  (Marine  Advisors,  Inc., 
1965**)  of  this  event  that  included  an  investigation  of  the  unpublished 
notes  of  the  late  Professor  G.  D.  Lauderback,  University  of  California, 

*  Multiply  feet  by  0.3048  to  convert  to  metres. 

**  Marine  Advisors,  Inc.  1965.  "An  Examination  of  the  Evidence  for  the 
Reported  Santa  Barbara  Coast  Tsunami  of  December  l8l2”  (unpublished 
report  submitted  to  the  Southern  California  Edison  Company). 


Berkeley  (cited  by  Wood  and  Heck  (1966)  as  a  basis  for  their  reported 
runup  heights)  has  shown  that  the  runup  heights  for  this  possible  tsunami 
were  not  nearly  as  large  as  those  reported.  Professor  Lauderback  took 
an  old  report  that  said  that  the  sea  flowed  1/2  mile  inland  and  trans¬ 
ferred  it  to  the  present  day  topographic  chart  in  the  vicinity  of  Gaviota 
Canyon.  This  yielded  a  50-ft  level.  However,  the  study  by  Marine  Ad¬ 
visors,  Inc.  (1965),  discovered  that  the  old  report  actually  referred  to 
an  estuary  near  Santa  Barbara  that  is  approximately  at  sea  level.  In 
addition,  they  learned  that  there  were  no  reported  drownings  even  though 
Indian  villages  were  on  the  beach  at  the  time  and  many  people  were  there. 
(The  mission  padres  kept  meticulous  records  of  Indian  births,  deaths,  and 
other  statistics).  Other  evidence  indicates  that  there  may  have  been  a 
runup  as  great  as  10  to  12  ft  at  Gaviota  and  a  smaller  runup  at  Santa 
Barbara  (Marine  Advisors,  Inc.,  see  footnote  on  page  7).  A  report  of  a 
tsunami  at  Santa  Cruz,  California,  in  18U0  also  has  been  shown  to  be 
erroneous . 

15.  The  largest  authenticated  locally  generated  tsunami  in 
southern  California  was  generated  by  the  1927  Point  Arguello  earthquake 
and  produced  runup  elevations  as  great  as  6  ft  in  the  immediate  vicinity. 
The  only  other  reports  of  locally  generated  tsunamis  in  southern  Califor¬ 
nia  were  for  possible  tsunamis  of  very  small  heights  in  Santa  Monica  Bay 
(Joy  1968). 

16.  Thus  relatively  small  locally  generated  tsunamis  have  been 
known  to  occur  in  southern  California,  but  there  are  no  reliable  reports 
of  major  locally  generated  tsunamis.  There  could  be  a  few  locations  in 
southern  California  (e.g.  Gaviota)  for  which  locally  generated  tsunamis 
are  more  significant  than  distantly  generated  tsunamis  because  the  eleva¬ 
tions  produced  by  distantly  generated  tsunamis  are  small.  However,  the 
frequency  of  occurrence  of  locally  generated  tsunamis  in  southern  Califor¬ 
nia  is  not  known  and  predictions  of  locally  generated  tsunami  elevations 
are  beyond  the  scope  of  this  report . 

17.  The  finite-difference  numerical  model  used  to  simulate  the 
generation  of  tsunamis  in  the  Alaskan  region  and  their  propagation  across 
the  deep  ocean  was  originally  developed  by  Hwang  et  al.  (1972),  and 
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described  in  detail  in  an  earlier  report  (Houston  and  Garcia  197^).  A 
numerical  model  developed  at  WES  (Garcia  1976)  that  used  a  more  efficient 
solution  scheme  (explicit)  than  the  Hwang  et  al.  (1972)  model  (implicit- 
explicit  scheme)  was  used  to  simulate  the  generation  of  tsunamis  in  the 
Peru-Chile  region  and  their  propagation  across  the  deep  ocean.  This 
model  was  needed  as  a  result  of  the  extremely  large  grids  (approximately 
100,000  grid  cells)  used  to  propagate  tsunamis  over  the  great  distance 
from  Peru-Chile  to  southern  California.  In  this  study,  the  model  employed 
1/5°  by  1/3°  spherical  coordinate  grids  to  solve  the  linearized  long-wave 
equations.  Such  coarse  grids  adequately  resolve  the  very  long  tsunami 
wavelengths  in  the  deep  ocean.  The  grids  covered  very  large  sections  of 
the  Pacific  Ocean  including  either  Alaska  and  the  west  coast  of  the  United 
States  or  South  America  and  the  west  coast  of  the  United  States.  The 
boundary  condition  on  the  solid  boundaries  (land)  of  the  grids  was  that 
the  component  of  the  velocity  normal  to  the  boundary  equaled  zero.  On 
open  boundaries  (ocean),  a  first-order  approximation  of  total  transmission 
was  made. 

18.  The  deep-ocean  finite-difference  model  solved  an  initial  value 
problem  starting  with  an  uplift  deformation  of  the  water  surface  identi¬ 
cal  with  the  major  features  of  the  permanent  deformation  (permanent  in 
the  sense  that  the  time  scale  associated  with  it  is  much  longer  than  the 
period  of  the  tsunami)  of  the  sea  floor  following  the  seismic  disturbance. 
The  transient  movements  within  the  time-history  of  the  ground  motion  were 
neglected  because  Hammack  (1972)  has  further  shown  that  the  initial 
deformation  of  the  water  surface  will  closely  approximate  major  features 
of  the  permanent  deformation  of  the  ocean  floor,  provided  these  features 
have  characteristic  lengths  that  are  at  least  four  times  as  great  as  the 
water  depth.  The  neglect  of  smaller  features  is  unimportant  because  such 
small-scale  details  produce  waves  that  are  negligible  in  the  far  field. 

Nearshore  Numerical  Model 


19.  The  nearshore  finite-difference  numerical  model  used  in  this 
study  is  described  in  detail  by  Butler  (1978)-  The  model  solves  the 
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classical  shallow- water  wave  equations  that  govern  tsunami  propagation 
in  the  nearshore  region.  A  highly  efficient  implicit  solution  scheme 
that  employs  a  centered,  alternating-direction  procedure  is  used. 

20.  The  nearshore  propagation  numerical  model  employed  in  this 
study  uses  a  smoothly  varying  grid  that  allows  cells  to  be  small  in  shal¬ 
low  water  and  larger  in  deeper  water.  A  piecewise  reversible  transforma¬ 
tion  (analogous  to  that  used  by  Wanstrath  (1976) )  is  used  independently 
in  the  x-  and  y-directions  to  map  the  grid  into  a  uniform  grid  used  in 
computational,  space.  The  coordinate  transformation  has  the  form 

.  r 

x  =  p  +  qot 

where  p  ,  q  ,  and  r  are  arbitrary  constants  and  the  transformation  is 
applied  piecewise  for  each  axis.  The  transformation  is  such  that  there 
is  a  regular  grid  in  computational  space  and  all  derivatives  in  computa¬ 
tional  space  are  centered.  Many  stability  problems  commonly  occurring 
in  variable  grid  schemes  are  eliminated  using  this  transformation  since 
the  real  space  grid  is  smoothly  varying  with  the  variation  and  its  first 
derivative  being  continuous.  The  variable  grid  allows  accurate  resolu¬ 
tion  of  the  tsunami  wave  form  as  it  enters  shallow  water  and  its  wave¬ 
length  decreases. 

21.  The  nearshore  numerical  model  used  as  input  the  time-history 
calculated  by  the  generation  and  deep-ocean  propagation  numerical  model. 

A  tsunami  was  generated  in  the  Aleutian-Alaskan  area  or  the  west  coast 
of  South  America  and  propagated  across  the  deep  ocean  to  a  500-m  depth 
off  the  west  coast  of  the  United  States.  Wave  forms  calculated  at  this 
depth  by  the  deep-ocean  numerical  model  were  recorded  all  along  the  west 
coast.  These  wave  forms  then  were  used  as  input  to  the  nearshore  numeri¬ 
cal  model,  which  propagated  the  tsunamis  from  the  500-m  depths  across  the 
continental  slope  and  shelf  to  shore. 

22.  The  nearshore  numerical  grid  covered  an  area  approximately 
280  miles  by  1^2  miles.  It  extended  from  Point  Conception,  California, 
to  Ensenada,  Mexico.  The  grid  was  oriented  approximately  parallel  to  the 
500-m  contour,  since  refraction  would  have  bent  the  wave  fronts  to  such 


an  orientation  before  they  reached  the  500-m  contour.  Grid  cells  in 
water  depths  less  than  500-m  were  slightly  less  than  1  mile  on  a  side. 

Verification 


23.  The  numerical  models  used  in  this  report  were  verified  by 
numerical  simulations  of  the  1964  Alaskan  tsunami.  This  was  the  only 
large  tsunami  for  which  reliable  information  exists  concerning  source 
characteristics.  The  initial  condition  used  in  the  generation  and  deep- 
ocean  propagation  numerical  model  was  that  the  uplift  of  the  ocean's 
surface  in  the  source  region  was  identical  with  the  permanent  deformation 
of  the  ocean  bottom  following  the  earthquake.  The  permanent  deformation 
of  the  ocean's  bottom  as  a  function  of  spatial  location  was  taken  from 
Plafker  (1964).  The  generation  and  deep-ocean  propagation  numerical 
model  was  used  to  propagate  the  1964  Alaskan  tsunami  to  a  water  depth  of 
approximately  500  m  off  the  coast  of  southern  California.  The  nearshore 
numerical  model  then  was  used  to  propagate  the  tsunami  to  shore. 

24.  There  were  seven  tide  gages  in  operation  along  the  coast  of 
southern  California  from  Point  Conception  to  San  Diego.  Figures  1-7  show 
comparisons  between  recordings  of  the  1964  Alaskan  tsunami  by  these  tide 
gages  and  wave  forms  calculated  using  the  generation  and  deep-ocean  propa¬ 
gation  and  the  nearshore  numerical  models  to  propagate  the  1964  tsunami 
from  Alaska  to  southern  California.  The  comparisons  extend  over  a  period 
of  1.5  to  2.0  hr.  This  time  period  was  selected  since  numerical  models 
that  generate  tsunamis  and  propagate  them  across  the  deep  ocean  to  distant 
locations  will  provide  good  values  only  for  the  first  1.5  to  2.0  hr 
(Houston  1978). 

25.  The  numerical  model  calculations  agree  remarkably  well  with 
the  tide  gage  recordings  shown  in  Figures  1-7,  considering  the  fact  that 
the  uplift  used  to  generate  the  1964  tsunami  is  only  approximately  known. 

In  addition,  tide  gages  do  not  record  tsunamis  without  .distortion  of 
amplitude,  period,  and  phase.  Even  so,  the  amplitudes,  periods,  and  phases 
of  the  calculated  and  recorded  wave  forms  shown  in  Figures  1-7  are  in  good 
general  agreement.  The  greatest  amplitude  difference  is  for  the  second 
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PART  III:  METHODOLOGY  FOR  ELEVATION  PREDICTIONS 


Tsunami  Occurrence  Probabilities 


26.  Historical  data  on  tsunami  generation  must  be  the  basis  for 
an  analysis  that  considers  the  probability  of  tsunami  generation  in  the 
two  tsunamigenic  areas  in  the  Pacific  Ocean  of  concern  to  southern 
California — the  Aleutian  and  Peru-Chile  Trench  areas.  A  satisfactory 
correlation  between  earthquake  magnitude  and  tsunami  intensity  has  never 
been  demonstrated.  Not  all  large  earthquakes  occurring  in  the  ocean 
even  generate  noticeable  tsunamis.  Furthermore,  earthquake  parameters 
of  importance  to  tsunami  generation,  such  as  focal  depth  and  vertical 
ground  motion,  have  only  been  measured  for  earthquakes  occurring  in 
recent  years.  Therefore,  data  on  earthquake  occurrence  cannot  be  used 
to  determine  occurrence  probabilities  of  tsunamis.  Historical  data  of 
tsunami  occurrence  in  generation  regions  must  be  used  to  determine 
these  probabilities. 

27.  In  South  America,  a  wealth  of  information  exists  concerning 
tsunami  occurrence.  Reliable  data  (grouped  in  intensity  increments  of 
0.5)  exist  for  tsunamis  with  intensity  greater  than  or  equal  to  0  for 
a  169-yr  period  and  greater  than  or  equal  to  2.5  for  a  Ul7-yr  period. 

The  intensity  scale  used  was  a  modification  (by  S.  L.  Soloviev  of  the 
Soviet  Union)  of  the  standard  Imamura-Iida  tsunami  intensity  scale. 
Intensity  is  defined  as 


1  ■  los2  (V5*^) 

This  definition  in  terms  of  an  average  runup  (in  metres)  over  a  coast 
instead  of  a  maximum  runup  elevation  at  a  single  location  (used  for  the 
standard  Imamura-Iida  scale)  tends  to  eliminate  a  spurious  intensity 
magnitude  caused  by  often  observed  anomalous  responses  (due,  for  example, 
to  local  resonances)  of  single  isolated  locations. 

28.  Using  the  most  recent  and  complete  catalog  of  tsunami  occur¬ 
rence  in  the  Pacific  Ocean  (Soloviev  and  Go  1969),  a  relation  between 
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tsunami  intensity  and  frequency  of  occurrence  was  determined  for  the 
tsunami-generating  trench  running  the  length  of  the  Peru-Chile  coast. 
Tsunamis  with  intensity  greater  than  or  equal  to  0  were  considered.  It 
was  assumed  that  the  logarithm  of  the  tsunami  frequency  of  occurrence 
was  linearly  related  to  the  tsunami  intensity.  Earthquake  magnitude 
and  frequency  of  occurrence  have  been  similarly  related  by  Gutenberg 
and  Richter  (1965)  and  used  extensively  in  earthquake  predictions. 
Soloviev  (1970)  has  shown  a  similar  relation  between  tsunami  intensity 
and  frequency  of  generation  for  moderate  to  large  tsunamis  throughout 
the  Pacific  Ocean.  Furthermore,  Wiegel  (1965)  found  the  same  type  of 
relation  between  tsunami  occurrence  and  runup  levels  for  historical 
tsunamis  at  Hilo,  Hawaii;  San  Francisco,  California;  and  Crescent  City, 
California;  and  Adams  (1970)  for  tsunamis  at  Kahuku  Point,  Oahu.  A 
recent  study  by  Rascon  and  Villarreal  (1975)  revealed  this  same  relation 
for  historical  tsunamis  on  the  west  coast  of  Mexico  (data  from  1732)  and 
on  the  Pacific  west  coast  of  North  America,  excluding  Mexico. 

29.  Letting  n(i)  equal  the  probability  of  a  tsunami  with  an 
intensity  i  being  generated  during  any  given  year  and  using  statistics 
for  the  entire  trench  along  the  Peru-Chile  coast,  a  least-squares  anal¬ 
ysis  resulted  in  the  following  expression: 

n(i)  =  0.07l+e~0,63i  (l) 

30.  In  using  statistics  for  the  entire  trench  area  along  the 
Peru-Chile  coast,  it  was  assumed  that  the  probability  of  tsunami  occur¬ 
rence  was  uniform  along  the  trench.  This  is  a  standard  assumption  for 
earthquake  frequency  analysis  (Gutenberg  and  Richter  1965 )•  The  tec¬ 
tonic  justification  of  this  assumption  lies  in  the  fact  that  a  single 
sialic  block  or  plate  of  the  earth’s  crust  or  lithosphere  is  dipping 
into  the  Peru-Chile  Trench  (Wilson  1 969)*  It  can  be  reasonably  expected 
that  the  movement  of  this  single  plate  is  similar  along  its  entire 
length. 

31.  In  the  Aleutian  Trench  area,  only  large  tsunamis  occurring  in 
relatively  recent  years  (since  1788)  have  been  recorded  due  to  the 
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isolation  of  the  area.  Assuming  an  exponential  coefficient  of  -0.71  for 
this  trench  area  (determined  by  Soloviev  1970)  as  a  mean  value  for  areas 
of  the  Pacific  with  the  most  data  on  tsunamis  and  using  only  the  re¬ 
liable  data  for  large  tsunamis  (intensity  greater  than  or  equal  to  3.5)» 
the  following  relation  was  determined  by  a  least-squares  analysis: 

n(i)  =  0.113e"0*71i  (2) 

Again,  the  probability  of  tsunami  occurrence  was  assumed  to  be  uniform 
along  the  trench. 


Use  of  Numerical  Models 


32.  To  relate  the  probability  distribution  of  tsunami  intensities 
to  source  characteristics,  it  was  assumed  that  the  ratio  of  the  source 
uplift  heights  producing  two  tsunamis  of  different  intensities  (as  de¬ 
fined  in  the  previous  section)  is  equal  to  the  ratio  of  the  average 
runup  heights  produced  on  the  coasts  near  these  tsunami  sources.  This 
ratio  is  equal  to  2(i.^  -  1^)  for  two  tsunamis  with  intensities  i^ 
and  i2  . 

33-  If  H  is  the  wave  height  emitted  in  the  direction  parallel 

8, 

to  the  major  axis  of  length  a  by  a  tsunami  source  with  an  elliptical 
shape  (large  tsunamis  have  historically  had  elliptically  shaped  uplifts) 
and  is  the  wave  height  emitted  in  the  direction  parallel  to  the 

minor  generation  axis  of  length  b  ,  then  experimental  research  of 
tsunami  generation  has  shown  that  =  a/b  (Hatori  1963).  For  a 

large  tsunami,  can  be  larger  than  H&  by  a  factor  of  as  much  ais 

5  or  6.  Thus,  the  orientation  of  the  tsunami  source  relative  to  the 
area  where  elevations  are  to  be  determined  is  very  important;  that  is, 
the  elevation  at  a  distant  site  due  to  the  generation  of  a  tsunami  at 
one  location  along  a  trench  cannot  be  considered  as  being  representative 
of  all  possible  placements  of  the  tsunami  source  in  the  entire  trench 
region.  Hence,  the  Aleutian  and  Peru-Chile  Trenches  had  to  be  segmented 
and  elevations  along  the  coast  of  southern  California  determined  for 
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tsunami  sources  located  at  the  center  of  each  of  the  segments. 

3**.  The  spatial  size  of  a  tsunami  source  was  standardized  because 
there  is  not  an  apparent  correlation  between  tsunami  intensity  and 
spatial  size  of  a  tsunami  source.  For  example,  the  19^6  Aleutian  tsu¬ 
nami  had  an  uplift  region  of  very  small  spatial  extent,  whereas  the  1957 
Aleutian  tsunami  had  an  uplift  region  that  covered  perhaps  the  greatest 
spatial  extent  of  any  known  earthquake  (Kelleher  et  al.  197*+);  yet  the 
19^6  tsunami  had  the  greater  tsunami  intensity,  producing,  in  general, 
greater  runup  elevations  in  the  near  and  distant  regions. 

35-  The  standard  source  employed  (described  in  Houston  and  Garcia 
197*0,  represents  a  large  tsunami  with  intensity  *+  on  the  modified 
Imamuia-Iida  scale.  Certainly,  tsunamis  of  low  intensity  may  have 
smaller  spatial  extents;  however,  large  tsunamis  pose  the  greatest 
threat  to  a  distant  area  such  as  southern  California.  These  large  tsu¬ 
namis  can  be  expected  to  have  similar  spatial  extents,  with  any  spatial 
differences  being  unimportant  in  the  far  field  compared  with  the  effects 
of  source  orientation  and  vertical  uplift. 

36.  Figures  8  and  9  show  the  Aleutian  Trench  divided  into  12  seg¬ 
ments  and  the  Peru-Chile  Trench  into  3  segments.  The  segments  in  the 
Aleutian  Trench  were  approximately  one-quarter  the  length  of  the  major 
axis  of  the  standard  source,  whereas  the  segments  in  the  Peru-Chile 
Trench  were  approximately  tne  length  of  the  major  axis  of  the  standard 
source.  The  standard  source  was  centered  in  each  segment  such  that  the 
major  axis  of  the  source  was  parallel  to  the  trench  axis.  Uplift  re¬ 
gions  historically  have  had  such  an  orientation  relative  to  trench  sys¬ 
tems.  The  Aleutian  Trench  was  segmented  much  finer  than  the  Peru-Chile 
Trench  because  the  Aleutian  Trench  is  oriented  relative  to  southern 
California  such  that  elevations  produced  in  southern  California  are  very 
sensitive  (Houston  et  al.  1975)  to  the  location  of  a  source  along  the 
trench.  Uplifts  along  the  Peru-Chile  Trench  do  not  radiate  energy  di¬ 
rectly  toward  southern  California  regardless  of  their  position  along  the 
trench.  Therefore  elevations  in  southern  California  are  not  very  sensi¬ 
tive  (Garcia  1976)  to  source  location  within  these  sections. 
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Idealized  axis  of  Aleutian  Trench  showing  12  segments  and  perimeter  of  standard  uplift  area 
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37-  In  each  of  the  segments  of  the  Aleutian  and  Peru-Chile 
Trenches,  tsunamis  with  intensities  from  2.5  to  4.5  in  increments  of 
0.5  were  generated  and  propagated  across  the  deep  ocean  using  the  deep- 
ocean  numerical  model  discussed  in  an  earlier  section.  Tsunamis  with 
intensities  less  than  2.5  are  too  small  to  produce  significant  eleva¬ 
tions  in  southern  California.  An  upper  limit  of  4.5  was  chosen  because 
the  largest  tsunami  intensity  ever  reported  was  less  than  4.5  (Soloviev 
and  Go  1969).  Gutenburg  and  Richter  (1965)  indicate  that  there  is  an 
upper  limit  to  the  strain  which  can  be  supported  by  rock  before  frac¬ 
ture.  Thus,  earthquakes  only  reach  certain  maximum  magnitudes  and  tsu¬ 
namis  can  be  expected  to  have  analogous  upper  limits  of  intensity. 

Perkins  (1972)  and  McGarr  (1916)  have  demonstrated  that  future  earth¬ 
quakes  cannot  have  seismic  moments  (measure  of  earthquake  magnitude  for 
large  earthquakes)  much  larger  than  those  of  earthquakes  that  occurred 
in  recorded  history. 

38.  The  wave  forms  propagated  to  southern  California  by  the  deep- 
ocean  propagation  numerical,  model  were  used  as  input  to  the  nearshore 
propagation  numerical  model.  Each  wave  form  was  propagated  from  a  water 
depth  of  500  m  to  shore  using  the  nearshore  model.  Thus,  at  each  grid 
location  on  the  shoreline  of  the  west  coast,  there  was  a  group  of  75  wave 
forms — 5  wave  forms  (for  intensities  from  2.5  to  4.5  in  increments  of 
0.5)  for  each  segment  of  the  Aleutian  and  Peru-Chile  Trenches.  Each  of 
these  wave  forms  had  an  associated  probability  equal  to  the  probability 
that  a  certain  intensity  tsunami  would  be  generated  in  a  particular  seg¬ 
ment  of  a  trench  region. 

Effect  of  Astronomical  Tides 


39.  The  maximum  "still  water"  elevation  produced  during  tsunami 
activity  is  the  result  of  a  superposition  of  tsunami  and  astronomical 
tide  wave  forms.  Therefore,  the  statistical  effect  of  astronomical 
tides  on  total  tsunami  runup  must  be  included  in  the  predictive  scheme 
presented  in  this  report.  Since  the  wave  forms  calculated  by  the 
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near-shore  model  did  not  have  a  simple  form  (e.g.,  sinusoidal),  the 
statistical  effect  of  the  astronomical  tide  on  tsunami  runup  had  to 
he  determined  through  a  numerical  approach. 

UO.  The  wave  forms  calculated  by  the  nearshore  numerical  model 
extended  over  a  period  of  time  of  approximately  2  hr.  One  to  three  wave 
crests  (the  largest  waves  in  the  tsunami)  arrived  during  this  time. 
Smaller  waves  arriving  at  later  times,  however,  have  often  persisted  for 
days  during  historical  tsunamis.  An  analysis  of  tide  gage  records  of 
the  i960  and  1961+  tsunamis  in  southern  California  indicated  that  these 
smaller  waves  have  amplitudes  on  the  average  of  50  percent  of  the  maxi¬ 
mum  wave  amplitude  of  the  tsunami;  therefore  a  sinusoidal  group  of  these 
small  waves  was  added  to  each  of  the  calculated  wave  forms  so  that  the 
total  wave  form  extended  over  a  12-hr  period.  These  smaller  waves  are 
important  for  locations  where  tsunami  waves  are  fairly  small  compared 
with  tidal  variations.  At  such  locations,  the  maximum  combined  tsunami 
and  astronomical  elevation  occurs  during  the  maximum  tidal  elevation. 

Ul.  A  computer  program  was  developed  to  predict  time-histories  of 
the  astronomical  tides  in  southern  California.  The  program  was  based 
upon  the  harmonic  analysis  methods  used  in  the  past  by  the  National 
Ocean  Survey  for  mechanical  tide-predicting  machines  (Schureman  19b8). 
Tidal  constants  available  from  the  National  Ocean  Survey  were  used  as 
input  to  the  computer  program.  A  year  of  tidal  elevations  was  then 
predicted  for  grid  locations  in  southern  California.  The  year  1961+  was 
selected  because  all  the  major  tidal  components  for  tides  in  southern 
California  had  a  node  factor  of  approximately  1.00  during  this  year, 
thus  making  it  an  average  year.  The  node  factor  is  associated  with  the 
revolution  of  the  moon’s  node  and  has  an  l8.6-yr  cycle.  Since  a  tsunami 
can  arrive  at  any  time  during  this  18.6-yr  period  (arrival  at  a  low  of 
the  node  factor  is  equally  likely  as  an  arrival  at  a  high),  the  statis¬ 
tical  effect  of  the  temporally  varying  node  factor  on  the  predicted 
runup  elevations  is  shown  to  be  very  small  in  Appendix  A  of  this  report. 

U3.  The  many  maximum  elevations  with  associated  probabilities 
were  used  to  determine  exceedance  frequency  distributions  of  combined 


tsunami  and  astronomical  tide  elevations.  The  maximum  elevations  were 
ordered  and  frequencies  summed,  starting  vith  the  largest  elevations, 
until  a  desired  frequency  was  obtained.  The  elevation  encountered  when 
the  summed  frequency  reached  a  desired  value  F  was  the  elevation  that 
is  equaled  or  exceeded  with  an  average  frequency  of  once  every  l/F  yr. 
Thus,  when  the  summed  frequencies  reached  the  value  0.01,  the  elevation 
associated  with  the  last  frequency  summed  was  the  100-yr  elevation. 


PART  IV:  RESULTS 


Explanation 

44.  Figures  10-23  present  predicted  100-  and  500-yr  elevations 

(in  feet)  in  southern  California  produced  by  distantly  generated  tsunamis. 
These  elevations  include  the  effects  of  the  astronomical  tide;  that  is, 
they  are  maximum  elevations  due  to  the  superposition  of  tsunami  and  tidal 
wave  forms  (see  PART  III).  The  lower  curves  in  Figures  10-23  represent 
the  100-yr  elevation  and  the  upper  curves ,  the  500-yr  elevation.  A 
100-yr  elevation  is  one  that  is  equaled  or  exceeded  with  an  average 
frequency  of  once  every  100  yr  (i.e.  100-yr  return  period);  a  500-yr 
elevation  has  a  corresponding  definition.  Elevations  in  this  report  are 
referenced  to  the  mean  sea  level  (msl)  datum. 

45.  The  elevations  presented  in  Figures  10-23  extend  along  more 
than  200  miles  of  coastline  in  southern  California  from  Santa  Barbara  to 
San  Diego.  Similar  elevations  north  and  west  of  Santa  Barbara  are  pre¬ 
sented  in  an  earlier  report  (Houston  and  Garcia  1978).  Calculations  were 
made  at  240  locations  (presented  in  Figures  24-31)  along  this  coastline. 
The  latitude  and  longitude  of  each  of  these  gage  locations  are  presented 
in  Table  1.  The  elevations  calculated  at  each  gage  location  are  connected 
by  straight  lines  in  Figures  10-23. 

46.  Plates  1-120  present  exceedance  distributions  at  all  240  gage 
locations.  The  y-axis  of  the  plots  presents  tsunami  frequency  in  events 
per  year.  Tsunami  frequency  is  the  inverse  of  the  tsunami  return  period. 
Thus  the  100-yr  tsunami  has  a  frequency  of  0.01  events  per  year  and  the 
500-yr  tsunami  a  frequency  of  0.002  events  per  year.  The  tsunami  fre¬ 
quency  is  presented  in  Plates  1-120  on  a  logarithmic  scale.  Each  tick 
mark  from  a  frequency  of  0.001  to  0.01  represents  an  increment  of  0.001. 
Each  tick  mark  from  a  frequency  of  0.01  to  0.1  represents  an  increment  of 
0.01.  Figure  32  illustrates  how  to  determine  100-  and  500-yr  elevations 
using  an  exceedance  frequency  distribution. 

47.  The  tsunami  elevations  presented  in  Figures  10-23  and 
Plates  1-120  are  elevations  at  the  shoreline.  As  discussed  by  Houston 
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Figure  15.  Gage  numbers  91-109 
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Figure  18.  Gage  numbers  1L5-163 
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Figure  19.  Gage  numbers  163-181 
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Figure  20.  Gage  numbers  181-199 
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Figure  29.  Gage  locations  158-193 


Figure  31.  Gage  locations  200-2140 
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EXCEEDANCE  FREQUENCY  DISTRIBUTION 

Figure  32.  Use  of  frequency  curves 


and  Garcia  (1974)  and  Houston  (1979) >  runup  elevations  should  he  equal 
to  elevations  at  the  shoreline  for  most  locations  in  southern  California. 
Thus  the  elevations  presented  in  this  report  also  can  be  considered  to 
be  runup  elevations.  However,  runup  elevations  may  not  equal  shoreline 
elevations  at  locations  where  dimes  prevent  flooding  except  through  in¬ 
lets  or  if  the  land  is  extremely  flat  and  inland  flooding  is  extensive. 

A  land-flooding  numerical  model  (e.g.  model  presented  by  Houston  and 
Butler  (1979))  can  be  used  to  determine  runup  for  such  cases. 

Discussion 


1+8.  The  exceedance  frequencies  presented  in  Plates  1-120  are  for 
combined  tsunami  and  astronomical  tide  elevations  given  the  occurrence 
of  a  noticeable  tsunami  in  southern  California.  It  is  not  true  that  in 
the  limit  of  a  tsunami  with  an  amplitude  of  zero,  the  frequency  analysis 
presented  determines  the  maximum  possible  elevation  for  the  astronomical 
tide  alone  over  a  100-yr  period.  It  is  quite  possible  that  the  100-yr 
combined  tsunami  and  astronomical  tide  elevation  is  less  than  the  maximum 
possible  astronomical  tide  elevation  over  a  100-yr  period  or  even  less 
than  the  maximum  astronomical  tide  elevation  that  occurs  during  a  typical 
year.  For  example,  the  maximum  combined  tsunami  and  astronomical  tide 
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elevation  during  the  last  100  yr  at  Rincon  Island,  California,  probably 
occurred  during  the  1964  Alaskan  tsunami  (the  I960  Chilean  tsunami  may 
have  been  somewhat  larger,  but  it  was  not  recorded  at  Rincon  Island). 

The  maximum  combined  elevation  was  approximately  3.5  ft  above  msl.  This 
elevation  is  greater  than  mean  higher  high  water  (2.8  ft)  but  less  than 
the  maximum  astronomical  tide  elevation  that  typically  occurs  during  the 
average  year  at  Rincon  Island.  Thus  the  fact  that  the  combined  tsunami 
and  astronomical  tide  elevations  presented  in  Plates  1-120  are  at  some 
locations  less  than  the  maximum  possible  or  even  yearly  tide  elevations 
that  may  occur  at  the  locations  is  a  property  of  the  natural  phenomenon 
(tsunami  heights  are  small  compared  with  the  tide,  and  arrival  of  a  tsu¬ 
nami  during  the  maximum  astronomical  tide  of  a  year  is  more  than  a  100-yr 
event)  and  not  of  the  methods  used  to  calculate  exceedance  frequencies. 

49*  The  logarithm  of  tsunami  frequency  is  plotted  versus  elevation 
in  Plates  1-120.  There  are  sections  of  the  exceedance  frequency  dis¬ 
tribution  curves  that  are  approximately  linear.  This  linear  relation¬ 
ship  between  the  logarithm  of  the  tsunami  exceedance  frequency  and  tsunami 
elevations  has  been  noted  by  many  observers  (see  PART  III). 

50.  The  linear  relationship  between  the  logarithm  of  tsunami 
frequency  and  tsunami  elevations  holds  only  over  a  certain  range  of 
tsunami  frequency.  Cox  (1964)*  has  shown  that  tsunami  frequencies  and 
elevations  are  linearly  related  for  frequencies  greater  than  approxi¬ 
mately  0.05  to  0.1.  Since  earthquakes  probably  only  reach  certain  maxi¬ 
mum  magnitudes,  tsunamis  can  be  expected  to  have  similar  upper  limits  to 
intensity.  Thus  the  exceedance  frequency  curves  must  eventually  approach 
a  maximum  elevation. 

Comparison  with  Predictions  Based  on  Local  Observation 

51.  Virtually  all  of  the  coast  of  southern  California  lacks  ob¬ 
servations  of  tsunami  activity.  Only  four  locations  (Alamitos  Bay,  Los 

*  Cox,  D.  C.  1964.  "Tsunami  Height-Frequency  Relationship  of  Hilo" 
(unpublished  paper),  Hawaii  Institute  of  Geophysics,  University  of 
Hawaii,  Honolulu,  Hawaii. 
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Angeles,  Santa  Monica,  and  San  Diego)  have  tide  gage  recordings  for 
both  major  tsunamis  in  i960  and  1964.  Since  the  tide  gage  at  Alamitos 
Bay  malfunctioned  during  the  i960  tsunami  and  there  is  a  gap  in  the 
recording  of  approximately  2-1/2  hr,  a  reasonable  determination  of  ex¬ 
ceedance  frequency  based  on  observed  data  cannot  be  made  at  Alamitos 
Bay.  The  tide  gage  at  Santa  Monica  also  malfunctioned  during  the  i960 
tsunami,  but  there  is  a  good  indication  of  a  lower  limit  for  the  maxi¬ 
mum  elevation. 

52.  Except  at  San  Diego  (where  a  tide  gage  has  operated  since  the 
midnineteenth  century),  recordings  of  tsunamis  in  southern  California 
were  not  made  prior  to  the  twentieth  century.  Thus  the  80-yr  period  of 
time  from  1900  through  1979  will  be  used  to  estimate  tsunami  exceedance 
frequencies  at  Los  Angeles,  Santa  Monica,  and  San  Diego  based  upon  local 
historical  data.  The  i960  and  1964  tsunamis  produced  the  largest  elevr- 
tions  at  all  three  of  these  locations  during  the  80-yr  period.  By  using 
data  for  these  two  tsunamis  and  assuming  that  the  logarithm  of  the 
tsunami  exceedance  frequency  is  linearly  related  to  the  maximum  combined 
tsunami  and  astronomical  tide  elevation  (see  PART  III),  the  100-  and 
500-yr  combined  tsunami  and  tide  elevation  can  be  estimated.  This  esti¬ 
mate  is  quite  approximate  since  the  time  period  is  relatively  short, 
there  are  many  locations  in  source  regions  where  tsunamis  can  be  gener¬ 
ated  but  were  not  during  the  time  period  (the  numerical  predictions  in¬ 
clude  this  effect  by  considering  tsunamis  generated  all  along  the  source 
regions),  and  the  effect  of  tsunami  arrival  during  various  tidal  stages 
is  not  accounted  for  in  the  historical  data. 

53-  Using  microfiche  of  tide  gage  recordings  (from  the  National 
Oceanic  and  Atmospheric  Administration)  of  the  i960  and  1964  tsunamis, 
the  i960  and  1964  tsunamis  were  found  to  have  produced  combined  tsunamic 
and  astronomical  tide  elevations  of  4.0  ft  and  2.8  ft  at  Los  Angeles, 
greater  than  4.6  ft  (tide  gage  went  off  scale)  and  3.0  ft  at  Santa  Monica, 
and  3.9  ft  and  3.1  ft  at  San  Diego.  The  tabulation  on  the  following  page 
compares  100-  and  500-yr  elevations  using  these  data  and  elevations  pre¬ 
sented  in  Plates  52,  38,  and  113  for  gage  locations  104  (Los  Angeles), 

75  (Santa  Monica),  and  225  (San  Diego)  respectively. 


Local  Historical 
Data  Predictions 
ft 


Predictions  Pre¬ 
sented  in  this 
Report ,  ft 


Location 

100  yr 

500  yr 

100  yr 

500  yr 

Los  Angeles 

4.4 

7.2 

5-2 

8.2 

Santa  Monica 

>5.2* 

>8.9* 

5.8 

10.2 

San  Diego 

4.2 

6.0 

4.6 

5.7 

*  Lower  limit 

since  tide 

gage  off 

scale  during  i960 

tsunami  at  Santa  Monica. 


PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


5*+.  The  deep-ocean  and  nearshore  numerical  models  used  in  this 
report  accurately  simulated  tsunami  generation,  deep-ocean  and  nearshore 
propagation,  and  interaction  with  coastlines.  The  nearshore  numerical 
model  was  superior  to  the  analytical  techniques  used  for  nearshore 
propagation  in  a  previous  study  (Houston  and  Garcia  197*0  because  it  was 
two-dimensional  and  handled  arbitrary  incident  wave  forms.  The  tech¬ 
niques  used  in  this  report  to  predict  tsunami  elevations  provided  good 
estimates  of  100-  and  500-yr  tsunami  elevations  for  southern  California. 
Based  upon  the  techniques  presented  in  this  report,  100-  and  500-yr 
predictions  at  Los  Angeles,  Santa  Monica,  and  San  Diego,  California, 
were  shown  to  agree  quite  well  with  predictions  based  upon  historical 
observations  of  tsunami  activity  at  these  three  locations.  These  were 
the  only  locations  in  southern  California  which  had  sufficient  local 
historical  data  to  allow  predictions  based  upon  local  observations. 

Recommendations 


55.  The  elevations  predicted  in  this  report  are  at  the  shoreline 
but  can  be  assumed  to  equal  runup  elevations  for  most  of  southern 
California.  There  are  locations  where  time-dependent  effects  (for 
example,  lack  of  sufficient  time  to  completely  flood  extensive  low- 
lying  or  estuarine  areas)  or  two-dimensional  effects  (for  example,  flow 
divergence  or  convergence)  cause  tsunami  runup  elevations  not  to  be 
equal  to  elevations  at  the  shoreline.  It  is  recommended  that  for  these 
areas  inundation  limits  be  determined  using  a  numerical  model  developed 
recently  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (Houston 
and  Butler  1979).  This  model  is  capable  of  handling  land  flooding  for 
bays,  harbors,  developed  areas  such  as  cities,  large  low-lying  areas, 
sand-dune  protected  areas,  and  other  areas  where  there  are  topographical, 
roughness,  or  coastline  variations. 
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Table  1  (Concluded) 
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APPENDIX  A:  EFFECT  OF  NODE  FACTOR  TEMPORAL  VARIATION 
ON  JOINT  PROBABILITIES 


1.  In  constructing  exceedance  frequency  distributions  for  ex¬ 
tremal  heights  resulting  from  Joint  occurrences  of  tsunamis  and  the 
astronomical  tide,  a  period  of  time  for  the  analysis  must  be  selected. 
However,  the  inclination  of  the  moon's  orbit  to  the  plane  of  the  earth's 
equator  (called  an  obliquity  factor)  varies  throughout  an  l8.6-yr 
cycle  because  of  the  revolution  of  the  moon's  node.  The  ratio  obtained 
by  dividing  the  true  obliquity  factor  by  its  mean  value  is  called  a  node 
factor  since  it  is  a  function  of  the  longitude  of  the  moon's  node.  The 
node  factor  thus  has  an  18.6-yr  cycle.  If  a  period  of  time  less  than 
18.6  yr  is  selected  for  an  analysis  of  joint  occurrences  of  tsunamis 
and  the  astronomical  tide,  the  influence  of  the  temporal  variation  of 
the  node  factor  is  neglected.  In  this  appendix  it  will  be  shown  that  if 
a  particular  year  (instead  of  18.6  yr)  is  selected  for  the  analysis, 

the  effect  of  the  temporal  variation  in  the  node  factor  on  the  exceedance 
frequency  distribution  is  negligible.  Thus  only  a  single  year  instead 
of  18.6  yr  needs  to  be  considered  in  the  analysis  and  this  results  in 
a  significant  reduction  of  required  computational  expense. 

2.  The  effect  of  a  temporal  variation  of  the  node  factor  on  ex¬ 
ternal  heights  resulting  from  joint  occurrences  of  tsunamis  and  the  as¬ 
tronomical  tide  can  be  estimated  using  an  analytical  solution  described 
by  Houston  and  Garcia  (19T^)»  Since  tsunamis  have  periods  much  less 
than  18.6  yr,  the  variation  in  the  node  factor  is  assumed  to  be 
negligible  during  a  tsunami  (in  fact,  the  node  factor  in  standard  tidal 
harmonic  analysis  is  assumed  to  be  constant  over  an  entire  year). 

Let 

z  =  the  water  elevation  at  any  time  above  local  mean  sea  level. 

P,p(z)  =  the  exceedance  frequency  distribution  for  the  water  level 
at  a  site  being  equal  to  or  exceeding  z  due  only  to  the 
maximum  wave  of  the  tsunami. 

P^(z)  =  the  exceedance  frequency  distribution  of  the  water  level  at 
a  site  being  equal  to  or  exceeding  z  due  only  to  the 
temporal  variation  of  the  node  factor  of  the  astronomical 
tide. 


P(z)  =  the  exceedance  frequency  distribution  of  the  water  level 
at  a  given  site  being  equal  to  or  exceeding  z  due  to 
both  the  maximum  wave  of  the  tsunami  and  the  node  factor 
of  the  astronomical  tide. 

3.  According  to  Chandrasekhar  (1943),  P(z)  can  be  calculated 


P(z)  =  J  fN(X)PT(z  -  A )dA 


where 


- 


and  fN(z)  is  the  frequency  density  for  the  astronomical  tide.  If 
PN(z)  is  approximated  by  a  Gaussian  distribution,  then 


t,{<)  =  -i- 
N 


where  the  variance,  o  ,  is  given  by 


-1  c2 


and  C  .  equals  the  height  variation  of  the  node  factor  associated  with 
m 

the  tidal  constituent. 

U.  It  is  well-known  that  the  logarithm  of  PT(z)  is  linearly 
related  to  tsunami  elevations  (see  PART  III).  Therefore,  let 


PT(z)  =  Ae" 


h 


A2 


Substituting  Equations  A3  and  A5  into  Equation  A1  yields 


-a 

P(z)  =  Ae 


£z-(oca2/2)J 


(A6) 


or 


P(z)  =  Ae  az 


1 


where 


(AT) 


(A8) 


5.  Thus,  the  net  effect  of  the  temporal  variation  in  the  node 
factor  is  to  produce  a  P(z)  identical  with  Pm(z)  except  for  a  shift 

2  2  1 

of  z  by  an  amount  acr  /2  .  The  variance  a  can  be  determined  by 
multiplying  the  nodal  factor  variation  for  each  tidal  constitutent 
(Schureman  I9I+8)  times  the  amplitude  of  the  tidal  constituents  (tidal 
constituents  available  from  the  National  Ocean  Survey).  This  product 
equals  in  Equation  A*+;  a  can  be  determined  at  a  location  using 

Plates  1-120  by  approximating  the  curves  with  straight  lines  and  deter¬ 
mining  the  slope  of  the  lines. 

6.  For  example,  gage  IOH  is  the  site  of  the  Los  Angeles  tide  gage. 

2 

From  Plate  52  a  is  found  to  be  equal  to  approximately  0.5*+ .  cj  equals 
0.17. 

2 

7.  Thus  ao  /2  £  0.05  ft  .  Hence,  the  exceedance  frequency  dis¬ 
tribution  based  upon  a  year  for  which  the  major  tidal  constituents  have 

an  average  value  (approximately  1.0)  should  be  shifted  relative  to  an 
l8.6-yr  period  by  approximately  0.06  ft  for  all  frequencies. 

8.  The  analytic  solution  presented  was  shown  to  be  correct  by  a 

direct  numerical  approach.  The  year  196*+  was  selected  as  the  year  for 

which  the  major  tidal  constituents  have  an  average  value.  The  l8.6-yr 
period  was  centered  on  196*+.  Tsunami  and  tidal  elevations  were  expressed 
•  ‘.he  nearest  0.1  ft.  It  was  found  that  the  exceedance  frequency  dis- 

’  .‘lens  for  the  1-yr  and  18.6-yr  period  had  the  same  50-yr,  200-yr, 
and  1000-yr  elevations.  The  100-yr  elevation  was  different  by 


A3 


V  > 

pT<  ) 


APPENDIX  B:  NOTATION 

Length  of  major  axis  of  elliptical  source,  miles 
Amplitude  constant 

Length  of  minor  axis  of  elliptical  source,  miles 
m  nodal  constituent 
Probability  density  function 
Tsunami  frequency 

Wave  height  in  direction  of  major  axis  of  ellipse,  ft 

Wave  height  in  direction  of  minor  axis  of  ellipse,  ft 

Average  runup  over  a  coast,  m 

Tsunami  intensity 

Tsunami  probability  function 

Constant 

Exceedance  frequency  distribution  for  combined  node  factor 
variation  and  tsunamis 

Exceedance  frequency  distribution  for  node  factor  variation 

Exceedance  frequency  distribution  for  tsuanmis 

Constant 

Constant 

Spatial  coordinate,  ft 
Spatial  coordinate,  ft 

Water-surface  elevation  above  mean  sea  level  (msl),  ft 
Exponential  constant 
Variable  of  integration 
Variance 
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